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PREFACE 

This is one of a continuing series of reports of the Ford Founda- 
tion sponsored Research Program in University Administration at the Uni- 
versity of California, Berkeley. The guiding purpose of this Progrcm is 
to undertake quantitative research which will assist university admini- 
strators and ether individuals seriously concerned with the management of 
university systems both to understand the basic functions of their com- 
plex systems and to utilize effectively the tools of modern management 
in the allocation of educational resources. 

This study investigates the resource allocation problem of faculty 
hiring and promotion patterns using the techniques of optimal control 
theory. If is both an extension and a synthesis of the conceptual ana- 
lysis of faculty structure introduced in earlier papers in this series* 
The principal characteristics of the faculty system considered are: 
(1) linear system propagation; (2) a convex preference function to rank the 
relative values of varying the states of the system; and (3) four state var- 
iables and four control variables including the stocks and flows of (a) full 
professors, (b) associate professors, (c) assistant professors, and (d) in- 
structors. The specific approach adopted for this investigation assumes 
that the promotion policies and attrition rates of faculty members are 
relatively fixed over the short run and the only variables left open to 
achieve a desired faculty structure are the institutional hiring policies. 
Under these conditic the optimal open loop faculty hiring paths are cai-« 
culated and their se^ditivity is investigated. Finally, this study inves- 
tigates and evaluates several solution procedures. 

This paper was presented at the Eleventh American Meeting of The In- 
ERJC stitute of Management Science, October 19-21, 1970. 



I. INTRODUCTION 

The analytical techniques of optimal resource allocation have been ap- 
plied to decision making in general and in institutions of higher education 
in particular."^ One area that is both financially dominant in the public 
sector and academically essential is the educational administrator's need 
to dec'-de on the most desirable pattern of faculty hiring and promotion. 

On a national basis, faculty salaries are the most expensive component of 

2 

an institution's budget, accounting for $5.5 billion in 1968-69, and the 
nu&iber of faculty positions is often fixed by law or slowly evolves under 
tenure restrictions. Furthermore, the quality of the faculty is essential 
because they set the tenor of the institution and participate in many of 
the academic and administrative decisions. 

Some aspects of this problem were discussed in a recen'c paper concer- 
ning equilibrium faculty patterns resulting from current in zutional ap- 

3 

pointment and promotion policies. Using the Berkeley campus of the Uni- 
versity of California as an example, it was shown that the logical extrapo- 
lation of the current data would result in no feasible equilibria for cur- 
4 

rent policies. More generally, present policies have resulted in facul- 
ties more heavily concentrated in tenure ranks than is often desired. Fur- 
thermore, under present policies this situation will tend to worsen over 
time because of the current youthful faculty age distribution reflecting 



For an extensive discussion of this research area, see Weather sby [16]. 

2 

Howard R. Bowen, "Financial Needs of the Campus," in Robert H, 
Connery [3]. 
3 

Oliver [ 9]. See also: Bartholomew [l], Halpern [ 4], Reisman [10], 
and Reisman and Taft [ll]. 
4 

In this context, feasibility was defined as those regions ^hich fulfilled 
^ certain tenure/non- tenure proportional constraints as well as the enforced dy- 

namics of the system under current appointment and promotion policies. 
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the rapid faculty expansion of the past decade. Therefore, one major pro- 
blem of the educational administrator Is to allocate his open faculty posi- 
tions wisely by choosing hiring and promotion policies which fulfill his 
long-term goals of faculty rank distribution while meeting his budgetary 
or other resource constraints. 

One approach to the analysis of this problem is for a campus or system 
level administrator to assume that in the short run the promotion policies 
and attrition rates of faculty members are both relatively fixed and beyond 
his lamedlate control. While it is true that over time an administrator 
can affect promotion and attrition rates by policy changes and financial in- 
centives, often the only control variables available to the campus decision 
maker to help him achieve a desired faculty structure in the short run are 
the institutional hiring policies. The scarce resource constraint can be 
Interpreted as either a limit on the total funds available for academic 
salaries or a limit on the total number of full time equivalent teaching 
positions allowed. 

While there are clear political and bureaucratic costs associated with 
exceeding an administrator's available resources, there is also concern as- 
sociated with the under-utill^atlon of a resource in the public sector. 
Unused re&ources often have high opportunity costs associated with them 
and public fiducial responsibility requires the maximum productive use of 
public funds. Furthermore, unmet public needs can foster political dis- 
content when funds are not used to their full availability. Finally, the 
future budget allocation of a scarce resource is frequently dependent upon 
the full use of that resource In the current period. Administrators are 
aware that under-utllization in one period may very well lead to a lesser 
budget for the following period* 

In addition to the resource constraint, the decision maket must also 



consider institutional preferences for a well-balanced faculty > not only in 
terms of subject field but also in terms of professorial rank. While the 
definition of this factor may vary widely, the balance of senior and junior 
faculty is never the lees important and relevant in the decision process. 

The decision system described in this study is a multi-stage decision 
system, in which the results of current decisions are perceived in subsequent 
years and current decisions are constrained by past hirii^g decisions. The 
degree to which the d*:*.cision maker fulfills his goals can be measured by a 
utility or scoring function which is defined for the particular decision 
maker in question. 

This paper is concerned with the application of control theory to the 
solution of .the optimal open-loop faculty hiring problem. Chapter II con- 
tains a description and estimation of the formal model of the faculty struc- 
ture. Chapter III contains a mathematical formulation of the analysis. 
Chapter IV contains discussion of the numei:ical results derived from the 
application of the model, and Chapter V summarizes the conclusions drawn 
from this study, and discusses future research. The computer program used 
for these solutions was written in a generalized format so that it could 
be easily applied by other researchers to similarly defined allocation pro- 
blems. Appendix B described how a copy of the user's manual or the pro- 
gram may be obtained frou: the authors. 
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II. THE CONCEPTUAL liODEL, DATA, AND ESTIMATION 
6H THfi SlSien DYNAMICS 

Conceptual Fonnulation 

For the purposes of this analysis, the variables which characterise 
faculty structure can be divided into (a) those variables which the deci- 
sion maker can directly control, called control variabilis and designated by 
the syiiibol u ; (b) those endogenous variables which cannot be controlled 
by the decision maker, called state variables and designated by the symbol 

X ; and (c) those exogenous variables impinging upon the system, which are 
designated by the symbol z . To avoid the possibility of the problem be- 
coming degenerate, we assume that neither x nor u is empty. ^ For the 
purposes of this analysis, chese variables are assumed to be related by 
the linear dynamic propagation equation: 

x(i + 1) - Fx(i) + Gu(i) + Hz(i) i « 0, 1, . i . , N - 1 (1) 

where the initial x(0) is given and where i is the planning period and 
N is the total number of planning periods considered. We eithei: kiiew or 
can measure the initial state of the system, x(0) . Also basic to the model 
is a preference function formulated by the decision maker to describe the 
relative values of the states of the system. The intensity of preference 
for a particular state can be written as a function of the current state var- 
iables, control variables and the decision period, and is denoted as 
V (x, u, i) . Bttcause w« ara eoncerned viteh an N period decision pro- 



id 
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If u is empty, che decision maker has no control over the faculty struc- 
ture and the decision problem is meaningless. If x is empty, the problem is 
directly analagous to the problem of consumer demand with the same solution. 
See Samuelson [12]. 
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blem, we must further define a scalar summary measure for the future stream 
of preferences,^ which we denote as J : 

N-1 

J -a^v(x(N), u*(N), n] + I a,v(x(l), u(l), l} . (2) 

1-0 ^ 

The N^*^ term Is separated from the first N-1 terms to reflect the 
truncation of an infinite sequence after N periods. The symbol de- 
notes a weighting or discount factor introduced to refle:t the time prefer- 
ence of the decision maker's utility. Using this notation, the multi-stage 
public sector resource allocation decision problem is to choose from the ad- 
missible set the control sequence ud), 1 ■ 0, 1, 2, . . . , N - 1 to 

N-1 

Max {j - a„V[x(N>, u*(N), N] + J] a V(x(i), u(i), il> (3) 
" 1-0 ^ 

subject to 

x(i + 1) - Fx(i) + Gu(i) + Hz(i) 1 « 0, ^ , . . . , N - 1 (4) 
x(0) given , 
(^and any budgetary or physical constraints. 

Application to the Faculty Structure Problem 

To apply this resource allocation decision model to the optimal faculty 
structure problem, we define the system variables as follows. The state and 
control vectors, x and u , refer to the four academic instructional ranks 
of all disciplines: (1) full professor, (2) associate professoi , (3) assis* 
tant professor, and (4) instructor, where x(i) are the number of faculty of 
ei^ch rank continuing at the end of academic year 1 , and u(i) are the num- 
ber of faculty of each rank hired at the beginning of year 1 . We may write 

6 

This assumes the intertemporal separability of preference to enable 
feasible assessment and tractable solution. 



the dynamics of the faculty structure as 
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In other words, the number of full professors in the system at the end of next 
year is a function of the number of full professors completing this year (per- 
sistence), the number of associate professors completing this year (promotion), 
and the number of full professors hired this year. We assume that (1) no fac- 
ulty member is promoted more than one rank in each year, (2) no faculty mem- 
ber is demoted to a lower rank, and (3) no newly hired faculty member is 
promoted within the first year of his contract. There have been very few 
exceptions to these assumptions in the recent history of the University of 
California. 

There are several logical restrictions thct should be imposed on the 
elemenus of the F and G matrices. Each year some people leave the sys- 
tem because of death, retirement, or resignation; therefore. 



0 <^ g^^ <^ 1.0 for i « 1, 2, 3, 4, and 



4 . 



\ £ . < 1.0 for j - 2, 3, 
i»l ^ 

The assumption of linear system equations allows one to use the tech- 
niques of multiple linear regression to estimate statistically the elements 
of the F and G matrices. Data were available for the total number of 
full time equivalent (FTE) faculty by rank for the academic years 1962-63, 
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to 1967-68 and a summary of the data is given In Table 1 and a correlation 

matrix is given in Table 2. The number of new hires for each rank and year 

were subtracted out for use as an independent variable in the regression 

equations. The results of the ordinary least squares multiple regression 

on data for total University of California are shown in Table 3. While 

2 

the coefficients of multiple determination (R ) were very high (.96 to .99), 
five of the eleven coefficients violated the logical sign and magnitude 
restrictions. . This is not surprising because the ordinary least squares 
regression did not recognize any inequality coefficient constraints. 

The estimation of transition probabilities has been discussed exten- 
sively in the literature and a number of techniaues have been developed^ 
to avoid the difficulties of ordinary least squares. Most of \ese tech- 
niques are formulated for equations in which the transition probability 
is tho dependent variable rather than the coefficient applied to an inde- 
pendent or predetermined variable. However, one approach that recognizes most 
of the logical coefficient restrictions is quadratic programming which 
minimizes [x(t + 1) - Fx(t) - Gu(t)]^lx(t + 1) - Fx(t) - Gu(t)] . This 
quadratic programming estimation technique was used on the same data and 
produced the coefficient estimates given in Table 4. We observe that all 
sign and magnitude restrictions are met by these estimates. 

It Is Interesting to compare the results derived by ordinary least 
squares with the estimates computed by quadratic programming. Table 5 shows 
that the sum of squared residuals and the standard error as a percent of the 
mean derived by quadratic programming estimation are very close to the cor- 
responding quantities derived by ordinary least squares, except in the case 

''see Lee, Judge and Cain [ 8] for a comparison of five alternative esti- 
mation techniques; Zellner and Lee [171 for a d^ -ussion of Logit, Frobit, 
Gompit, weighted least squares and joint estimation techniques; and Theil 
O and Rey [13] for a discussion of quadratic programming. 

ERIC 
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of instructors. 

Finally, a set of subjectively assessed coefficients was used cs a test 
for sensitivity and reasonableness. These values are given in Table 6 and 
indicate that annually 10% of the instructors are ^. acted to assistant pro- 
fessors, and 20% of both the assistant professors and the associate professors 
are promoted one rank. Furthermore, theso Judgmental values indicate that 
only 5% of the faculty Ittcvetthe system each year. This is greater than the 
recent attrition experience of about 2% at the University of California. 
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TABLE 3: TRANSITION MATRICES DERIVED FROM UNCONSTRAINED 
MULTIPLE REGRESSIONS 



Full 

Associate 
Assistant 
Instructor 



Full 
0.6854 
0.0 
0.0 
0.0 



F-Matrix 

Associate 
1.11 
-0.0689 
0.0 
0.0 



Assistant Instructor 



0.0 

0.1088 
0.A956 
0.0 



0.0 

0.0 
-0,5818 
-0.2955 



Full 

Associate 
Assistant 
Instructor 



Full 
0.9112 
0.0 
0.0 
0.0 



6 -Matrix 

Associate 
0.0 
2.050 
0.0 
0.0 



Assistant 
0.0 
0.0 

C.8138 
0.0 



Instructor 
0.0 
0.0 
0.0 

0.2413 



Full 



Z-Vector ' constants) 



Associate Assistant Instructor 



-413.6 



703.9 



275.2 



28.26 
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TABLE 4: TRANSITION MATRICES DERIVED FROM QUADRATIC 
PROGRAMMING ESTIMATION 



Fiai 

Associate 
Assistant 
Instructor 



Full 
0.7058 

0.0 

0.0 

0.0 



F-Matrix 

Associate Assistant Instructor 
0.5242 0.0 0.0 



0.9570 



0.0 



0.0 



0.03 



0.960 



0.0 



0.0 



0.450 



0.526 



6-Matrix 



Full Associate Assistant Instructor 



Full 

Associate 
Assistant 

Instructor 



1.000 

0.0 

0.0 

0.0 



0.0 

0.63 

0.0 

0.0 



0.0 
0.0 
0.23 
0.0 



0.0 
0.0 
0.0 

0.738 
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TABLE 5: COMPARISON OF EFFICIENCY OF ORDINARY 
LEAST SQUARES AND QUADRATIC PR0GRAM1ING ESTIMATION 



Full 

Ordinary Lease Squares 

Sum o£ Squared Residuals 4960 

Standard Error o£ Estimate 49.8 

Standard Error as Percentage 

o£ Mean 3.04 

Quadratic Programming 

Sum o£ Squared Residuals 5045 

Standard Error o£ Estimate 50.2 

Standard Error as Percentage 
o£ Mean 3.07 



Associate Assistant Instructor 

11590 19690 730 
76.1 99.2 15.6 

8.30 7.30 15.97 

11591 33541 1513 
76.1 129.5 22.5 

8.30 9.50 22.95 
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TABLE 6: JUD6MENTALLY ASSESSED TRANSITION MATRICES 



F-Matri)i 



Full 

Associate 
Asslbtant 
Instructor 



Full 

0.95 

0.0 

0.0 

0.0 



Associate Assistant Instructor 



0.20 



0.75 



0.0 



0.0 



0.0 



0.20 



0.75 



0.0 



0.0 



0.0 



0.10 



0.0 



Full 

Associate 
Assistant 
Instructor 



Full 

0.90 
0.0 
0.0 
0.0 



G-Matrix 



Associate Assistant Instructor 



0.0 
0.90 
0.0 
0.0 



0.0 
0.0 
0.90 
0.0 



0.0 
0.0 
0.0 
1.0 
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III. FORMULATION OF THE ANALYSIS 



Chapter II presented the general framework of public sector resource al- 
location decision analysis in terms of the decision maker's utility function 
V defined over the state and control variables. However, it is often more 
convenient to transform the analysis from one of maximizing utility to one of 
minlaizing loss, with less defined around the most desired targets. Near the 
targets , the loss structure is approximately quadratic independent of the 

8 

form of the utility function. 

Several criterion functions were used and most combined considerations of 
the relative composition of faculty by rank and either a monetary or absolute 
constraint on the total number of taculty positions. The relative composi- 
tion of the faculty was measured by the ratio of Zhe number in each rank to 
the number of full professors. The desired distribution is then a target or 
goal sought in some or all periods. The various criterion functions are sum- 
marized in Table 7. 



^If V(x, u) has a maximum at x*, u* (desired targets), then the second 
order Taylor series expansion about x*, u* is ^ 

V(x, u) - V(x*.u*) + l^x + l^u + l/2Ax^ l^x + l/2Au^ |^ Au + 



Near the optimum 



t a^v 

+ Ax'^^V^" + HOT . 

dxdu 



3V 
3x 



3u 



and 



x*,u* x*,a* 



-V(x*,u*) + 1/2[Ax'Au*^] 


V 1 
XX 'xu 




Ax 




V V 




Au 




ux uu 





+ HOT 



Therefore, the maximization of V(x, u) is equivalent to minimization of the 
quadratic term on the right hand side. 
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TABLE 7 

SUMMARY OF CRITERION FUNCTIONS USED IN THIS ANALYSIS 



Criterion Function 


Targets Met in 


Constrained by 


Weights Used 


1 


all periods 


budget 


various 


2 




budget 


various 


3 


all periods 


budget 


Kalman 


4 


last period 


budget 


various 


5 


last period 


budget 


Kalman 


6 


alirppclods 


positions 


various 


7 


last period 


positions 


various 



Deietljig Pj (1) as the average annual salary for faculty of rank j 
in period 1 , B(l) as the total academic salary budget for period 1 , r^ 
as the target ratios sought, and and 3 as weights Indicating relative 

loss, we may write the criterion function used as 



4 



V(l) 



1-2 J 



x^(i) 



. ., N - 1 



j-2 

for i - 0, 1, 



+ B 



I p (i)(x (i) + u (i)) 

J-1 



- B(i) 



(1) 



and 



V(N) 



4 

j-2 J 



x (N) 

[ 1 



(2) 



The weights describe the relative loss of deviating from the distri- 

butional targets of each rank. Several methods were used to derive these 
weights including (1) the average salary of each rank relative to the average 
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salary of full professors expressing the relative monetary loss of a unit 
deviation from the desired targets, (2) the Kalman variance in each period 
of each state variable relative to the variance in that period of the number 
of full professors expressing the relative uncertainty in each state var- 
iable, and (3) weights chosen roughly proportional to the first derivative of 
V(i) to insure rapid convergence. With tha exception of the Kalman variance 
weights, both types (1) and (3) were used in most numerical experiments. 

There are other important objectives served in the management of a fac- 
ulty structure. In some cases, the total number of faculty positions is fixed 
in each period and all composition adjustments must occur within these ceil- 
ings. This case can be accommodated within our framework by setting all the 
annual faculty salaries equal to 1.0 and the total academic salary budget 
equal to the chosen ceilings. Two other objectives not included in this ana- 
lysis are (1) the maintenance of a steady flow of promotions to avoid faculty 
ossification and discouragement and (2) the elimination of no longer produc- 
tive faculty through early retirement or some other means. These last two 
examples show that this analysis cannot solve all of the problems in the manage- 
ment of an academic faculty structure. On the contrary, this study illustrates 
that some faculty management decisions can be analyzed in a cogent and sophis- 
ticated manner. 
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Stochastic Considerations 

In addition to the deterministic model Just described, we may also 
introduce a stochastic element In the academic structural equations in the 
following fashion. We consider an additive random error in a linear dynamic 
system: 

x(i + 1) - Fx(i) + Gu(i) + H2(i) + e (i) (3) 
where e (i) is assumed to be independently normally distributed with zero 

9 

mean and varlance-covariance matrix Q , which is assumed constant over time. 
The decision maker can observe the state of his system at the beginning of per- 
iod i by the linear scheme 

7(i) - S(i) x(i) + v(i) (4) 

where S(i) is the current samplirg scheme matrix, v(i) is the sampling er- 
ror with variance-covariance matrix R(i) , and y(i) the observation vector. 
It may now be derived^^ through use of the Kalnan filter that, while the prior 



Estimates of the elements of the variance-covarif.nce matrix Q can be de 
rived from the results of either the linear regressions or the quadratic pro- 
gramming estimation performed to develop the F and G matrices. Where e^^. 
is the i^h residual and ij the mean residual of the J^^ prediction equa- 
tion (i.e., for state J - 1.2, NX), and N is the number of data points 
used in each calculation, we have that 

Var(e, J - — n z for J - 1, 2, . . ., NX 

ij N - i 

and N 



Cov(e^j, e^j^) 

for j - 1, 2 NX ; k - 1, 2 NX ; j / k . 

^^Bryson and Ho [2], Chapter 12, 
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variance-covariance matrix of x before measurement, M , is 

M(i + 1) - FP(i)F^ + Q , (5) 

the posterior variance-covariance matrix of x after measurement, P , is 

P(i) - M(i) - M(i)S^(i)[S(i)M(i)S^(i) + R(i)]"''"S(i)M(i) . (6) 

In this way, given M(0) , Q , H(i) , and R(i) for i - 0, 1, . . ., N , 
we may precompute the posterior variance-covariance matrix, P , before attack- 
ing the deterministic optimization problem. While this stochastic element is 
not considered in the optimization process p^r 86, its use in the formulation 
of the preference function can be quite relevant. In the special case of a 
quadratic criterion function linear system dynamics perturbed by additive 
Gaussian noise, a "certainty equivalence" principle is applicable which al- 
lows one to separate the estimation and optimizaton procedures without af- 
fecting the final optimal solution.^"'' This is due to the fact that for the 
quadratic preference function, the expected value of the preference func- 
tion, V(i) , can be separated into its mean and variance. This is also 
true of negative-exponential and linear preference functions. 
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IV. NUMERICAL RESULTS 

This chapter is devoted to presenting numerical examples of the 
results derived from the Implementation of the model discussed In the 
previous chapters. All of the examples described here were designed to In- 
vestigate the behavior of the model under various conditions. Basic 
to this Investigation were the following characteristics: 

1. The criterion functions used were those described In Section 
III. 

2. All the faculty ratios are In proportion to the number of full 
professors. 

3. Five or ten planning periods were used. 

4. The data employed were for the total University of California. 

The Initial state variable assignments were adapted from the 

12 

actual FTE faculty appointment schedule for the 1^67^68 aca- 
demic year: full professors (1807), associate professors (822), 

13 

assistant professors (1189), and Instructors (13). 

5. No discounting was used on the preference function (i.e., all 

- 1.0). 

6. All prices and the total number of positions are assumed to 
Increase 5.0% per ye&r while the budget constraint Increases 
at the historical 12.0% per year. 

7. The control assignments of the Initial Iteration were all zeroes. 



8. For practical purposes, the computations were terminated when 
neither t\ie control set nor the criterion function value show 
any pronounced variation* This was usually 50-100 Iterations 



12 

University of California, 1967-68 Statistical Summary - Students and 
Staff [14]. 
13 

The unusually low figure for Instructors was due to the one-year con- 
tract for faculty of this rank. Practically all Instructors are new hires. 
14 

Q The Inaccuracies caused by any pre«-mature terminations were thought 

gPJ^(]; to be of minimal Importance for this particular study due to the cancella- 

— " tlons of error uoon comparison with other similarly handled runs. 
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In addition to examining several criterion functions, a Kunber of 
variations in the basic formulation were investigated by changing the 
relative size of the parameters of the criterion function. Three specific 
variations were in the: (1) price and budget vectors, (2) relative p-nalty 
weights, and (3) faculty ratios. In addition to these preference function 
variations, three different Sets of F and G transition matrices were 
used as were several different solution algorithms including both first 
and second order methods. Examples of these variations along with their 
derived results will now be described. 
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Analysis with Constrained Least Square Transition Matrices and Second 
Order Optimization Technique 

Budget Constrained 

The first investigation focused on criterion functions 1 and 4 which 
represent an institutional decision maker striving to achieve faculty 
distributional targets in each period or only in the final period while 
meeting the budget constraints. In this case, five planning periods were 
used. The results of thie analysis are given in Table 8 and Figures 1 
and 2. 

In all the analyses performed, the distribution of faculty chosen 
as the desired target was either the distribution included in the Regents 
Budget of the University of California, the existing pattern in 1967-68, 
or an arbitrary target. (See below) 

Table 8 



Source 



Distributional P^atios 



Actual 1967-68 
Budget 1967-68 
Arbitrary 



Assoc. Prof. 
Full Prof." 

.459 

.544 

.900 



Asst. Prof. 
Full Prof. 

.666 

1.192 

1.500 



Instr. 



Full Prof, 
.007 
.200 
.050 



Table 9 shows that even undsr optimal control the arbitrary faculty 
distributional targets were unachievable in any period because of resource 
and system constraints and there was very little difference in performance 

O 

ERIC 



between the two criterion functions. Furthermore, as in virtually all 
of the cases investigated, the budget constraint was very closely approxi- 
mated. 

Figure 1 shows the values of the control variables that are optimal 
for each of the two criterion functions under open loop control. While 
both hiring patterns average roughly 300 new faculty f each rank in each 
year, the patterns have striking differences. The Mrge number of new 
assistant professors hired under the first criterion function is necessary 
to maintain the requisite balance in the face of the severe estimated 
first year attrition of new assistant professors (see Table 4). Meanwhile, 
the fourth criterion function includes faculty distribution in the last 
period only and, consequently, the last decisions diverge considerably 
from the previous pattern in an effort to achieve the desired targets. 

From Figure 2 we observe that the total number of faculty of each 
rank also display quite different growth patterns. Under the first cri- 
terion function, the number of full professors increases very slowly at 
first and more rapidly at the end of the planning period, while under the 
fourth criterion furction the number of full professors follows exactly 
the opposite growth pattern. The number of instructurs is anothar extreme 
case which varies from a smooth increase (»/l) to an initial jump of 700 
followed by a gradual decline of 200 positions (if4X, It is interesting to note 
that with the same resources available and striving for the same targets 
an institution managed by criterion function 1 acquires more of the tenured 
ranks while an institution managed by criterion function 4 acquires propor- 
tionately more junior ranks. 
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Position Constrained 

The second series of analyses focuses on a similar set of decision 
problems with total position constraints instead of total budget con- 
straints. The faculty distributional targets in this case are the bud- 
geted ratios of 1967-68. Once again, the decision maker could seek to 
achieve these targets in every decision period (#6) or only in the final 
decision period (#7), The results of this analysis are given in Table 
10 and Figures 3 and 4. 

The evidence of Table 10 indicates that the budgeted faculty dis- 
tributional targets are also infeasible with the system description and 
position constraints used in this analysis- While the numbers of as- 
sociate professors exceed Its target, both assistant pi^ofessors and 
instructors have significantly lower ratios. • Meanwhile, both formu- 
lations closely adhere to the forecasted budgets. 

The initial presence of full professors in excess of the proportion 
indicated by the budgeted targets leads the analysis to recommend firing 
several hundred full professors as shown in Figure 3. While this is in- 
stitutionally infeasible and non-negativity conptralnts could be imposed 
on the new hires, these results are included to show the logical conse- 
quences of the budgeted faculty distributional ratios. The average new 
hires of about 100 per rank per year is significantly less than the ap- 
proximately 300 new hires of each rank in each year found under the bud- 
get constraint because the number of positions increases at an assumed 
5.0% per year as opposed to an assumed net 7,0% per year budget increase. 
The patterns of total faculty at each position are shown in Figure 4 
which again reveals abrupt alterations in the last decision period for 
the case of distributional targets in the last period only. 
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Use of Varianc s as Penalty Weights 

The weights used in the previous criterion functions were chosen either 
to balance the convergence of the solution algorithm or to reflect the 
relative cost of a one unit deviation from the targets. Another approach 
is to incorporate in the criterion function the current uncertainty in the 
estimate of the future magnitude of the state variables. In other words, 
a decision maker may choose to exert stricter control over the most uncer- 
tain components of his system. As discussed in Section III, the propa- 
gation of variance in a linear dynamic system is independent of the control 
chosen and, therefore, may be computed in advance and used as a set of 
fixed weights in the optimization. Furthermore, the relative magnitudes 
of the variance of the state variables will differ if the decision maker 
chooses to sample in the future. The results calculated using equation 
(III-6) for variance propagation are given in Table 11 for the case in 
which faculty distributional targets are sought in every decision period 
anJ in Table l2 for the case of only final period target ratios. The 
corresponding hiring and total faculty sequences are shown in Figures 5, 
6, 7, and 8. 

While the previous example penalized deviations from target ratios 
equally, the relative variances differ by about 50-100 to 1 with much less 
weight given the instructor target. As a consequence, the optimal number 
of newly hired instructors goes negative at some point in 3 out of A of 
the cases shown in Figures 5 and 7. This suggests chat an educatior^l 
decision maker would readily eliminate instructors to maintain balance in 
his other ranks - fe result observed in practice (see Table 8). Other- 
wise, the hiring patterns for the top three ranks with the variance weightings 
closely resemble the previous arbitrarily weighted results shown in Figure 1. 
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As before, the terainal tarcet only case exhibits more extreme behavior 
in the last decision period. 
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TABLE 13: SUMMARY OF THE EFFECTS OF CHANGING THE PLANNING HORIZON 

A decision maker with a preference function that seeks to achieve 
faculty ratios in last period only. 

assoc/fun asst/fun inst/full budget 

weights ; 1000. 0 500.0 250.0 0.0000005 



Period Salaries G^nerited; 10 Xbi Salfcrlet Budgeted Salarii 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 



88446.0625 
99073.875 
111155.875 
124507.4375 
139624.5000 
155578.1875 
175045.6250 
195762.8750 
218870.4750 
244766.5625 



88516.00 
99138.00 
111034.00 
124358.00 
139281.00 
155995.00 
174714.00 
195680.00 
219161.00 
245461.00 



Ratios 



assoc/fuH 



asst/fjll 



inst/full 



Targets 0.459 
Generated (10 Yr.) 0.519 
Generated' ( 5 Yr.) 0.475 



0.666 
0.617 
0.661 



es Generated; S Yr . 

88361.875 

99161.0625 
111467.500 
124746.125 
139421.4375 



0.007 
0.101 
0.064 



Preference Function Values 

10 Year Horizon : 

1) 0.002446 2) 0.002056 3) 0.007427 4) 0.01117 5) 0.05900 

6) 0.08687 7) 0.05499 8) 0.003434 9) 0.04223 10) 0.2411 11) 7.005 

Criterion Function Value : 1.251 

5 Year Horizon : 

1) 0.01188 2) 0.0002659 3) 0.09396 4) 0.07532 5) 0.009861 
6) 1.060 

Criterion Function Value; 1.251 
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Analysis of Sensitivity to Alternative Transition Matrices • First 
Order Optimization Technique 

Basic to the propagation of the state of the system over time is 
the structure and composition of the F and G transition matrices. 
The resulting effect of these matrices on the optimum allocation pattern 
is vividly shown in Figure 11, These curves show the results of using 
two different sets of transition matrices on otherwise identically defined 
allocation problems. Example (a) uses matrices derived from uncon- 
strained multiple regressions on actual data (see Table 3). Example (b) 
uses matrices assessed by subjective reasoning (see Table 6), 

The total faculty curves for example (a) show some system instability 
as opposed to the smooth growth pattern characteristic of example (b). 
This is due to the multiplicative effect of the coefficients present in the 
transition matrices estimated by regression analysis which contained pro- 
blems of realism, as earlier discussed. These results should be compared 
with Figure 2 (#1) which used the constrained least squares estimate. 

This comparison accentuates the importance of developing valid and 
reasonable coefficients to describe the propagation of the system, 
without which om cannot hope to produce reasonable results. It appears 
that constrained least squares regression is the best approach for determining 
these coefficients provided enough valid data are available. However, 
care must be taken in usifig the results of the regressions to evaluate 
their accuracy and reasonableness. 
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Variation of Weight on Resource Constraint with First Order Solution 
Technique 

As described in Section III, the preference function can be divided 
into two parts: (1) the faculty distribution constraint and (2) the total 
resource constraint. Each of these terms is weighted by coefficients 
which determine how strictly the constraint is to be enforced. A large 
weighting coefficient will increase the penalty resulting from the 
deviation; similarly, a small weighting coefficient will decrease the 
penalty. This section investigates the effects of the weighting coeffi- 
cient, 3 , on the resource constraint. 

Again we compare the results of decision problems which are identical 
in every respect except in the single characteristic of the budget weighting 
coefficient, 3 • These runs were made using the ordinary least squares 
regression results for transition matrices. Figure 12 shows the results 
of a moderate variation in 6 • The resource constraint is still operative 
when 3 is decreased from 0.01 to 0.001, although it is obviously not 
as effective. Notice that there is little change in the degree to which 
the system continues to fulfill the faculty ratio constratnt. 

In Figure 13, we see a much greater effect resulting from a variation 
in B . In thi s case, reducing 3 from 0.001 to 0.0001 makes the 
resource constraint completely inoperative. There is now a negative cor- 
relation between the budgeted total faculty and the calculated faculty. 
Furthermore, we note that there is now a marked increase in the system's 
ability to fulfill the faculty ratio constraints. 

In other words, the relative size of 3 regulates the balance between 
the resource constraint and the faculty ratio constraints. Consequently, 
3 must be chosen carefully to simulate correctly the goals of the 
decision maker. 
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Variation in Method for Calculating the Step-size in First Order Solution 
Procedure 

When investigating the convergence properties of the solution to 
this formulation, ; two different methods for calculating the incremental 
step-size, e , were investigated: (1) a single step-size, and (2) a 
multiple step-size. With the single- method, the same step-size is applied 
to the entire control set through the relationship 

new old - ou 

In this way, the size of the increment is proportional to the associated 
gradient. 

The multi-e method uses an individually computed step-size, ^^(i) > 
for each control eJ aent and is applied using the relationship. 

Under the multi-^: method, adaptive increments are computed by bisecfing 
the previous iteration's step-size when the associated gradient changes 
sign. 

Although the single-e method was the moat straightforward and most 
easily implemented of the two procedures, its convergence peoperties often 
proved to be less satisfactory than the convarg«nce of the maltl^ method. In 
alaost all cases » the speed of convetgeoce dropped off considerably as the opti- 
mum point was approached. This was caused by a large decrease in the size 



E V(x,u,t) + X^(t i- l)tFx(t) + Gu(t)]. 
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of the gradient which then resulted in too small an incremental step- 
size. Moreover, attempts to correct this by starting with a large initial 
e-value resulted in divergence of the solution algorithm. These problems 
with the single-e method led to the development of the multi-e method. 
The result was a much quicker overall convergence. 

An example of the difference in convergence properties for the two 
methods is shown in Figure 14. There is a break-even point between the 
speed of convergence for the two methods; consequently, a preferable mixed 
strategy would be to change from the single-e method to the multi-e method 
when the slopes of the two curves are equal. However, there is no method 
available at this time for determining where this point occurs. 
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Comparison of Three Solution Algorithms 

Finally, we present a comparison of the first-order gradient and 
multi-e methods and the second order gradient method. Table 14 compares 
similar solutions to one of the decision formulations and shows that the 
second order method is much better at meeting the chosen targets than 
either of the other methods, while maintaining approximately a 1% 
tolerance on the budget target. For this reason, the second order metho 
was primarily used on this analysis. 
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' TABLE 14: SUMMARY OF SOLUTION ALGORITHMS 

r 

MuUi-e Gradient Second Order 

Budget Weight .0001 .0000001 .000001 

Faculty Weights 100,50.25 1000,800,500 1000,800,500 

F and G matrices UNCONSTBAINED MULTIPLE REGRESSION 



After 30 Iterations with initial controls • 0 



Ratios 


Targets 


Achieved 


Achieved 


Achieved 


assoc/full 


0.9 


0.363 


0.582 


0.897 


asst/£ull 


1.5 


0.245 


0.426 


1.480 


inst/full 


0.05 


0.024 


0.016 


0.035 



Salaries Budget 

1 88516. 

2 99138. 

3 111034. 

4 124358. 

5 139281. 



% Diff. 

88507. .01 

99149. .01 

111039. .01 

124373. .01 

139276. .01 



% Diff . 

59102. 32.23 

57702. 41.8 

5/397. 48.31 

56401. 54.65 

55845. 59.91 



% Diff . 

89568. 1.19 

98171. .08 

110989. .05 

119774. 3.69 

139181. .08 



J -value 



107.7 



2728. 



25.32 
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V. CONCLUSIONS AND FUTURE RESEARCH 

This analysis has described the mathematical structure of an academic 
faculty by a linear dynamic model whose parameters were estimated from ac- 
tual data by two different techniques. The new faculty hiring decision pro- 
blem was formulated as an optimal control problem defined in terms of desired 
faculty distributional targets. Finally, either total budget or total posi- 
tion constraints were imposed upon the system through a penalty function. 
This multi-stage optimization problem subject to dynamic constraints was 
then solved by three alternate numerical methods^ 

While many interesting and illuminating -'esults were obtained and dis- 
cussed in this paper, a number of unresolved problems remain. The proper 
estimation of meaningful transition matrices is difficult without a long his- 
tory of individual faculty flows. The transition probabilities used in this 
analysis were based on aggregate data available for an entire institution 
and consequently obscured individual faculty time paths through the rank 
structure. If this kind of model is to be used for institutional decision 
making, a more complete and disaggregated data base would be needed. 

The criterion functions used in this analysis focused exclusively on the 
number and mix of faculty independent of student enrollments, available 
facilities, support staffs, or other relevant institutional parameters. This 
deliberate abstraction enabled a first step of demonstrable computational 
feasibility; however, the scope of the model should be expanded in future 
formulations to include at least the student body, the physical plant and 
the financial parameters of the system. 

Uncertainty is explicitly included in tha system dynamics in the form of 
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a random additive error term. This presumes that the transition probability 
matrices are known and fixed, which certainly is not the case. A full ana- 
lysis of uncertainty would include random transition matrices as well as ran- 
dom errors and would recognize the improved information available through 
sampling in the future - information which would provide new, updated esti- 
mates of transition probabilities. 

The conceptual formulation used in this analysis is most appropriate 
for a campus or system executive officer and provides little guidance to 
the ch-^irman of a department. Presumably, the department chairman can keep 
most of the relevant parameters in mind while he considers possible new fac- 
ulty members. However, a college chancellor or president must balance many 
competing pressures calling for more faculty and would probably desire a 
more disaggregate model than the one discussed in this paper. It would be 
feasible to extend the current model to the discipline or departmental level 
and to recognize additional categories of faculty including visiting and ir- 
regular faculty and teaching assistants. This would necessitate a larger 
computer program and a calculation time longer than the current 1-2 minutes 
on an IBM 360/ fi5; one hundred organizational units or groupings and ten cat- 
egories of faculty could easily be accomadacii* 

This analysis is only a beginning. It shows that a decision and con- 
trol formulation, which is more comprehensive and informative than simula- 
tion analysis, is economically feasible and potentially very useful in the 
analysis of faculty hiring decisions. The problems and difficulties of 
this approach deserve the attention of quantitative analysts. The promise 
and potential of this approach deserve the serious consideration of academic 
administrators. 
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APPENDIX A: SOLUTION METHOD 

The solution of this faculty allocation problem may be approached in 
a number of ways, one of which Is the use of the techniques of modern control 
theory. Since the mathematical background of these methods is adequately de- 
scribed elsewhere,^ only a brief discussion of the specific technique used 
will be presented here. 

We begin the solution method by adjoining the linear propogation con- 
straint to the criterion function with the U8« -6fi a Lugiange multiplier, 
thereby reducing the problem to an unconstrained optimization problem. 
Letting 

V(i) « V(x(i), u(i). i] (1) 

and 

f(i) - Fx(i) + Gu(i) + H2(i) (2) 

we may write the augmented criterion function, J , as 

N-1 

J - V(N) + I {V(i) + A^i + l)(f(i) - x(i + 1)1} • (3) 
i-0 



We now define the Hamiltonian, H(i) , as 
_H(i) « V(i) + A^(i + l)f (i) 



"''Bryson, A. Jr., and Yu-Chi Ho, Applied Optimal Control: Optimlza 
tion. Estimation and Control , (Waltham, Massachusetts; Blaisdell Publishing 
Co.), 1969. 

2 

For convex criterion functionsvlth I linear constraints the stationary 
solution of the adjoineu criterion function is the same as the stationary 
solution of the original criterion function. In other cases, the proper- 
ties of the criterion function and constraint set need to be investigated 
carefully. 



(4) 
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and may subsequently write the augmented criterion function as 

N-1 

J - V(N) + y [H(i) - x'(i + l)x(i + 1)1 . (5) 

i-0 

By changing the Index of suomation, we obtain 

N-1 

J - V(N) - X^N)x(N) + [ [H(i) - X^(i)x(i)] +H(0) . (6) 

1-0 

For a stationary point, we require that the total differential of J 

3 

be equal to zero along the entire sequence of allocations. Therefore, 



du(0) . (7) 



Since the initial state of the system, x(0) , is assumed given and fixed » 
dx(0) « 0 . In order to force dJ to zero we require that 

X^i) "H^iy i - 0, 1, 2 N-1 (8) 



subject to the boundary condition that 

av(N) 

3x(N) 



and also that 

0 i - 0. 1. . . N . 1 . (10) 

This provides us with the following first-order conditions necessary for an 
optimal allocation pattern: 



3 

The requisite second order conditions necessary for a maximum or a 
minimum are given in Bryson and Ho, op. cit., Chapter 6. 
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(1) x(i + 1) = Fx(i) + Gu(i) + Hz(i) i = 0, 1, . . . , N - 1 (11) 
x(0) given; 



(2) X^i) - i - 0, 1 N - 1 



and 



(3) 0 = 1^ i - 0, 1 N - 1 

It should be noted that nowhere in the previous discussion are there 
any specific limitations on the format of the preference function, 
V[x, u, i] , other than certain convexity requirements^ and the necessary 
inclusion of at least two non-empty variables sets: u(i) and x(i) . In 
other words, the fornjulation of the preference function is left completely 
up to the user. 



See footnote 2. 
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APPENDIX B: OPCON: A USER'S MANUAL 

A copy of both the uscr'i manual and the program are available upon re- 
quest and may be obtained by writing 

Dr. George B. Weathersby 

Ford Foundation Research Program 

Office of the Vice President - 

Planning and Analysis 
247 University Hall 
University of California 
Berkeley, California 94720. 

A nominal fee will be charged for duplicating and loading the program. The 
manual is available at no charge. Raquests ror further information should 
be sent to the same address. 
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